Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  594 1803 (N31E57 L=123)   884    108        1 ( 569/1773- 626/1833) -50.2 E60-E52 N27-N33
      2 1265 1684 (N40E04 L=175)   586     72        0 (1232/1650-1297/1718) -26.6 E07-E01 N38-N42
      3  635 1651 (N20E43 L=137)   367     60        1 ( 605/1637- 667/1670)  11.3 E45-E40 N19-N21
      4  702 1533 (N14E33 L=147)   587    165        1 ( 678/1471- 734/1607) -45.0 E37-E29 N10-N18
      5  366 1009 (S23E42 L=138)   414    117        1 ( 324/ 976- 415/1032)  33.7 E44-E37 S25-S21
      6 1460 1010 (N04W25 L=205)   643    195        2 (1376/ 986-1551/1032)  18.4 W19-W31 N01-N06
      7 1565  891 (N00W35 L=215)   294     59        2 (1534/ 868-1586/ 922)  56.3 W33-W36 S01-N02
      8  701  838 (S26E15 L=165)   517    190        0 ( 606/ 807- 781/ 873)  35.0 E20-E10 S29-S21
      9 1569  744 (S07W40 L=220)  1356    244        2 (1514/ 634-1611/ 839) -77.0 W36-W43 S14-S01
     10 1326  696 (S17W27 L=206)  1650    358        2 (1155/ 642-1494/ 773)  35.8 W16-W34 S22-S08
     11 1010  459 (S40W18 L=197)  1838    340        4 ( 851/ 427-1180/ 494)  16.3 W05-W29 S44-S35
     12 1514  398 (S27W54 L=233)  1750    389        2 (1350/ 367-1652/ 428)  29.7 W42-W63 S32-S20
  

    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