Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  385 1414 (N29E44 L= 17)  3920    456        8 ( 270/1246- 488/1609) -33.1 E76-E32 N17-N45
      2 1545 1483 (N22W41 L=101)   754    181        2 (1452/1456-1623/1500) -13.0 W34-W47 N19-N24
      3  167 1208 (N18E65 L=356)  2865    263        7 (  95/1089- 244/1288)  32.9 E77-E55 N11-N22
      4 1457  966 (S08W26 L= 87)   277     88        1 (1434/ 925-1473/1012)  68.2 W24-W27 S10-S05
      5  427  703 (S15E44 L= 17)   959    258        1 ( 355/ 602- 491/ 799)  34.5 E52-E36 S20-S09
      6  211  737 (S10E64 L=356)   415     64        0 ( 198/ 700- 229/ 773)  38.7 E67-E60 S12-S08
      7 1260  738 (S20W12 L= 73)   212     87        0 (1227/ 715-1311/ 760)  21.8 W10-W15 S21-S19
      8 1080  665 (S23W00 L= 60)   756    232        2 ( 986/ 635-1179/ 729) -32.5 E05-W06 S25-S18
      9  236  576 (S20E71 L=350)  2414    145        6 ( 186/ 505- 279/ 628) -45.0 E79-E65 S25-S17
     10 1462  461 (S42W29 L= 89)   915    127        2 (1414/ 429-1511/ 504) -11.3 W24-W34 S43-S38
     11 1273  417 (S43W13 L= 73)   265     45        1 (1242/ 393-1296/ 442) -26.6 W10-W14 S44-S40
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y